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^3 ' A new phase of a 116 Cd double j3 decay experiment is in progress in the 

< 



Solotvina Underground Laboratory. Four enriched 116 CdW04 scintillators with 



' total mass of 339 g are used in a set up, whose active shield is made of 15 

.— 

^ ; 2.5-3.2 MeV is 0.03 counts/yr-kg-keV. The half-life for 2v20 decay of 116 Cd is 

measured as T 1/2 (2u) = 2.6 ± 0.1(stat)l^(syst) xlO 19 yr. The T 1/2 limits for 

O ' 

1 neutrinoless 2/5 decay of 116 Cd are set at T 1/2 > 0.7(2.5) x 10 23 yr at 90%(68%) 

0) ' 

C.L. for transition to ground state of 116 Sn, while for decays to the first 2f and 
second 0^ excited levels of 116 Sn at T 1/2 > 1.3(4.8) x 10 22 yr and > 0.7(2.4) x 10 22 
yr with 90%(68%) C.L., respectively. For 0^2/3 decay with emission of one or 

><! ' 

two Majorons, the limits are T 1/2 (0z/Ml) > 3.7(5. 8)xl0 21 yr and T 1/2 (0i^M2) 



> 5.9(9.4) xlO 20 yr at 90%(68%) C.L. Restrictions on the value of the neu- 
trino mass, right-handed admixtures in the weak interaction, and the neutrino- 
Majoron coupling constant are derived as: m v < 2.6(1.4) eV, r/ < 3.9xl0~ 8 , 

A < 3.4xl0" 6 , and g M < 12(9.5) xlO" 5 at 90%(68%) C.L., respectively. 
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I. INTRODUCTION 



Neutrinoless (Ou) double /3 decay is forbidden in the Standard Model (SM) since it violates 
lepton number (L) conservation. However many extensions of the SM incorporate L violating 
interactions and thus could lead to the Oz/2/3 decay |],||]. Currently, besides conventional neu- 
trino (u) exchange mechanism, there are many other possibilities to trigger this process ||. In 
that sense neutrinoless 2/3 decay has a great conceptual importance due to the strong Schechter- 
Valle theorem || obtained in a gauge theory of the weak interaction that a non-vanishing 0^2/3 
decay rate requires neutrino to be massive Majorana particle, independent of which mecha- 
nism induces it. Therefore, at present 0v2f3 decay is considered as a very powerful test of new 
physical effects beyond the SM, and even the absence of this process - at the present level of 
sensitivity - would help to restrict or narrow this wide choice of theoretical models. At the 
same time Oz/2/5 decay is very important in the light of the measured deficit of the atmospheric 
neutrinos flux and the result of the LSND accelerator experiment which could be 

explained by means of the neutrino oscillations requiring in turn nonzero neutrino masses (m v ). 
However oscillation experiments are sensitive to neutrino mass difference, while measured (V2/? 
decay rate can give the absolute value of the m v scale, and hence provide a crucial test of m v 
models. 

Despite of numerous attempts to observe 0z/2/3 decay from 1948 up to present days this 
process still remains unobserved. The highest T"i/ 2 (0z/) limits were set in direct experiments 
with several nuclides: T 1/2 > 10 22 yr for 82 Se 100 Mo §, 116 Cd §; T 1/2 > 10 23 yr for 130 Te 
g and 136 Xe ]TT[; and T 1/2 > 10 25 yr for 76 Ge [0,0]. 

With the aim to enlarge the number of 2/3 decay candidate nuclides studied at a sensitivity 
comparable with that for 76 Ge and 136 Xe (neutrino mass limit of 0.5-2 eV), cadmium tungstate 
crystal scintillators, enriched in 116 Cd to 83%, were developed and exploited in 116 Cd research 
H4]||. The measurements were carried out in the Solotvina Underground Laboratory in a salt 



mine 430 m underground (~1000 m w. e.) |L5|. In the first phase of the experiment only 
one 116 CdWO"4 crystal (15.2 cm 3 ) was used. It was viewed by a photomultiplier tube (PMT) 
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through a light-guide 51 cm long and placed inside a plastic scintillator (038x115 cm) which 
served as veto counter. A passive shield of high purity (HP) copper (5 cm), lead (23 cm) and 
polyethylene (16 cm) surrounded the plastic counter. The background rate in the energy range 
2.7-2.9 MeV (Q 2/3 =2805 keV 0) was equal ^0.6 counts/yr-kg-keV. With 19175 h statistics the 



half-life limit for 0z/2/5 decay of 116 Cd was set as T 1/2 (0v) > 3.2xl0 22 yr (90% C.L.) f|, which 
corresponds to the restriction on the neutrino mass m u < 3.9 eV |J7]. Limits on 0^2/3 decay 



with emission of one (Ml) or two (M2) Majorons were obtained too: 7i/ 2 (0z/Ml) > 1.2xl0 21 
yr and T 1/2 (0z/M2) > 2.6xl0 20 yr (90% C.L.) p|. 



In the present paper new and advanced results of 116 Cd research obtained with the help of 
an upgraded apparatus are described. 

II. NEW SET-UP WITH FOUR 116 CDW0 4 DETECTORS 
A. Set-up and measurements 

In order to enhance the sensitivity of the 116 Cd experiment, the following improvements 
were scheduled: increase of the number of 116 Cd nuclei, reduction of the background and 
improvement of the data taking and processing |J. With this aim the new set-up with four 
enriched 116 CdW04 crystals (total mass 339 g) has been mounted in the Solotvina Laboratory 
in August 1998. All materials used in the installation were previously tested and selected for 
low radioactive impurities in order to reduce their contributions to background. 

In the new apparatus, a scheme of which is shown in fig. [I], four enriched crystals are viewed 
by the PMT (EMI9390) through one light-guide 10 cm in diameter and 55 cm long, which 
is composed of two glued parts: quartz 25 cm long and plastic scintillator (Bicron BC-412) 
30 cm long. The 116 CdW04 crystals are surrounded by an active shield made of 15 natural 
CdW04 scintillators of large volume (~200 cm 3 each) with total mass of 20.6 kg. Due to 
the high purity of the CdW04 crystals |19] and their large density (~8 g/cm 3 ) this active 



shield reduces background effectively. The veto crystals are viewed - by a low background 
PMT (017 cm) - through an active plastic light-guide (017x49 cm). In turn the whole array 
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of CdW04 counters is placed inside an additional active shield made of polystyrene-based 
plastic scintillator with dimensions 40x40x95 cm. Thus, together with both active light-guides 
(connected with enriched and natural crystals on opposite sides), a complete 47r active shield 
of the main 116 CdW04 detectors is provided. 

The outer passive shield consists of HP copper (thickness 3-6 cm), lead (22.5-30 cm) and 
polyethylene (16 cm). Two plastic scintillators (120x130x3 cm) are installed above the pas- 
sive shield to provide a cosmic muons veto. Because air in the Solotvina Laboratory can be 
contaminated by radon (at the level < 30 Bq/m 3 ) the set-up is isolated carefully against air 
penetration. All cavities inside the shield are filled up by pieces of plexiglass, and HP Cu shield 
is sealed with the help of silicon glue and enclosed inside a tight mylar envelope. 

The new event-by-event data acquisition is based on two IBM personal computers (PC) and 
a CAMAC crate with electronic units. For each event the following information is stored on 
the hard disc of the first computer: the amplitude (energy) of a signal, its arrival time and the 
following additional tags: the coincidence between different detectors; the signal of radio-noise 
detection system; triggers for light emitting diode (LED) and pulse shape digitizer. The second 
computer records the pulse shape of the 116 CdW04 scintillators in the energy range 0.25-5 
MeV. This complementary system is developed on the basis of a fast 12 bit ADC (Analog 
Devices AD9022) and is connected with computer by parallel digital I/O board (PC-DIO-24 
from National Instruments) [p2]| . Two PC-DIO-24 boards are used to link both computers and 
establish - with the help of proper software - a one-to-one correspondence between the pulse 
shape data recorded by the second computer and the information stored in the first PC. 

The energy scale and resolution of the main detector - four enriched crystals taken as a 
whole - were determined in the measurements with different sources ( 22 Na, 40 K, 60 Co, 137 Cs, 
207 Bi, 226 Ra, 232 Th and 241 Am). The energy dependence of the resolution can be expressed (for 
the energy above 50 keV) as FWHM(keV) =V-226 + 1Q.QE + 6.42 x 10^ 3 £ 2 , where energy 
E is in keV. For instance, the resolution (FWHM) was equal to 14.5% at 1064 keV and 11% 
at 2615 keV. The full energy peaks are well fitted in the energy region 0.06-2.6 MeV by a 
Gaussian function with typical value x 2 — 0.8-1.9. Moreover, the calibration spectra of the 
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232 Th source were simulated with the help of GEANT3.21 package |23] and event generator 
DECAY4 |24| (the last defines initial kinematics of the events). The simulated 232 Th spectra 
are in good agreement with the measured ones confirming the assumption of a Gaussian peak 
shape. In particular for the 2615 keV peak of 208 T1 - which is close to the 2/5 decay energy of 
116 Cd - the value of x 2 =1-3. 

Also, the relative light yield for for a particles as compared with that for electrons (a//3 ra- 
tio) and energy resolution were measured with a sources and corrected by using time-amplitude 
analysis (see below) as following: a/ (3 = 0.15(1) +7xl(T 6 £ a and FWHM a (keV) = 0.053£ Q 
(E a is in keV). The routine calibration is carried out weekly with a 207 Bi source (570, 1064 and 
1770 keV) and once per two weeks with 232 Th (2615 keV). The dead time of the spectrometer 
and data acquisition is monitored permanently with the help of an LED optically connected 
with the main PMT. The actual dead time value is ~4.2 % (~3% is owing to random coincidence 
between the main and shield detectors; ~1.2% is caused by miscounts of the data acquisition). 

The background spectrum measured during 4629 h in the new installation with four 
116 CdW0 4 crystals is given in fig. ^ where the old data obtained with one 116 CdW0 4 crystal 
of 121 g are also shown for comparison. As it is visible from this figure, the background is 
decreased in the whole energy range, except for the (3 spectrum of 113 Cd (Q@ = 316 keV), 
whose abundance in 116 CdW04 crystals is ~2% |14|]. In the energy region 2.5-3.2 MeV - where 



the peak of 0z/2/5 decay of 116 Cd is expected - the background rate is reduced to a value of 0.03 
counts/yr-kg-keV (only 4 events in the energy window 2.5-3.2 MeV were detected during 4629 
h), twenty times lower than in the previous set-up. It is achieved, first, due to improvement of 
passive and active shield, and secondly, as a result of data processing advance (time-amplitude 
and pulse-shape analysis), which are described below. 

B. Time-amplitude analysis of the data 

The energy and arrival time of each event can be used for analysis and selection of some 
decay chains in 232 Th, 235 U and 238 U families (see f. e. ref. pT|JT4| ). As an example (important 
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in the following for the background rejection in the energy range of Oz/2/3 decay), we consider 
here in detail the time-amplitude analysis of the following sequence of a decays from 232 Th 
family: 220 Rn (Q a = 6.40 MeV, T 1/2 = 55.6 s) -> 216 Po (Q a = 6.91 MeV, T 1/2 = 0.145 s) -> 
212 Pb. Because the energy of a particles from 220 Rn decay corresponds to ~1.2 MeV in 
scale of 116 CdW0 4 detector, the events in the energy region 0.7-1.8 MeV were used as triggers. 
Then all events (within 0.9-1.9 MeV) following the triggers in the time interval 10-1000 ms 
(containing 94.5% of 216 Po decays) were selected. The spectra of the 220 Rn and 216 Po a decays 
obtained in this way from data - as well as the distribution of the time intervals between 
the first and second events - are presented in fig. |3|. It is evident from this figure that the 
selected spectra and time distribution are in an excellent agreement with those expected from 
a particles of 220 Rn and 216 Po. Using these results and taking into account the efficiency of the 
time-amplitude analysis and the number of accidental coincidences (3 pairs from 218 selected), 
the determined activity of 228 Th ( 232 Th family) inside the 116 CdW04 crystals is as low as 38(3) 
//Bq/kg. 

The same technique was applied to the sequence of a decays from the 235 U family: 223 Ra 
(Q a = 5.98 MeV, T 1/2 = 11.44 d) -> 219 Rn (Q a = 6.95 MeV, T 1/2 = 3.96 s) -> 215 Po (Q a = 7.53 
MeV, Ti /2 = 1.78 ms) -> 211 Pb. For the fast couple ( 219 Rn -> 215 Po) all events within 0.8-1.8 
MeV were used as triggers, while a time interval 1-10 ms (65.7% of 215 Po decays) and an energy 
window 0.9-2.0 MeV were set for the second events. The obtained a peaks correspond to an 
activity of 5.5(14) /iBq/kg for the 227 Ac impurity in the crystals. 

As regard the 226 Ra chain ( 238 U family) the following sequence of j3 and a decays was an- 
alyzed: 214 Bi (Q p = 3.27 MeV, T 1/2 = 19.9 m) -> 214 Po (Q a = 7.83 MeV, T 1/2 = 164.3 fj,s) 
— > 210 Pb. For the first and second events the energy threshold was equal 0.1 MeV, and a time 
interval of 100-1000 /is (64.1% of 214 Po decays) was used. While the obtained spectrum of 
the first pulses agrees with the model of the (3 decay of 214 Bi, and the distribution of the time 
intervals between the first and second events can be fitted by an exponent with 7\/ 2 = 140l 2 g 
//s (in reasonable agreement with the 214 Po half-life value), the spectrum of the second events is 
continuous, contrary to the anticipated a peak of 214 Po. Probably, part of this continuous dis- 
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tribution can be explained by 226 Ra contamination of the materials neighboring the 116 CdW04 
crystals (optical grease, teflon, Mylar, radon in air), while another part is caused by 226 Ra 
decays in the crystals. Under such an assumption activity limits for 226 Ra contaminations are 
derived as <0.13(3) Bq/kg for optical grease, <8 mBq/kg for teflon, <1.8 /iBq/dm 2 for Mylar, 
and < 5 //Bq/kg for 116 CdW0 4 , whose values do not contradict bounds obtained earlier [|TS|j . 
To prove these assumptions, the events belonging to the 214 Bi — > 214 Po — > 210 Pb chain were 
independently searched for in the time window of 5-88 fis (28.9% of 214 Po decays) with the 
help of pulse shape analysis (see below). For both events the energy threshold was ~0.3 MeV. 
The result obtained ( 226 Ra activity in the 116 CdW04 crystals < 14 /iBq/kg) is similar to that 
of the time-amplitude analysis. 

Finally, all couples of events found for 232 Th, 235 U and 238 U families as described above were 
eliminated from the measured data. 



C. Pulse-shape discrimination 

The pulse shape of the 116 CdW04 scintillators in the energy region of 0.25-5 MeV is digitized 
by a 12 bit ADC and stored in 2048 channels with 50 ns channel's width. Due to different shapes 
of scintillation signal for various kinds of sources^ (a particles, protons, 7 quanta and cosmic 
muons were investigated), the pulse-shape (PS) discrimination method based on the optimal 
digital filter |20| was developed and clear discrimination between 7 rays (electrons) and a 



particles was achieved p2f . 

The pulse shapes of enriched crystals were measured for a particles with an 241 Am source 
and for 7 rays with 60 Co, 137 Cs, 207 Bi and 232 Th sources in the special calibration runsQ. To 



1 It is known, that scintillation efficiency and pulse shape of inorganic crystals depend on the local 



density of the energy released, hence allowing to identify the incoming radiation (see e. g. ref. [ 20 , 21j] ) . 

2 Because 7 rays interact with matter by mean of the energy transfer to electrons, it was assumed 
that pulse shapes for electrons and 7-s are the same. This statement was proved in the measurement 



7 



provide an analytic description of the a or 7 signals f a (t) and / 7 (t) the pulse shape resulting 
from the average of a large number of individual events has been fitted with the sum of three (for 
a particles) or two (for 7-s ) exponents, giving the reference pulse shapes f a (t) and fJt) (see 
for more details ref. p2|). In the data processing the digital filter is applied to each experimental 
signal f(t) with aim to obtain the numerical characteristic of its shape (shape indicator, SI) 
defined as: SI = J2k /(^fe) ' Pifk ~ to), where the sum is over all time bins (from k = 1 to k = 
2048), f(tk) is the digitized amplitude of a given signal (normalized to its area) at the time 
The weight function P(t k — t a ) is determined as P(t) = {/ ' a (t) — f 7 (t)}/{f a (t) + / 7 (£)}, and t a 
is the time origin of the signal. The measured with sources SI distributions are well described 
by a Gaussian functions, whose mean values and standard deviations a a and cr 7 have a slight 
energy dependence^. For 0.9 MeV 7 quanta S7 7 = 18 ± 3, while for 4.8 MeV a particles SI a 
= 29.0 ± 3.6. It allows us to determine the efficiency of the PS event selection for the different 
chosen intervals of SI values (±cr, ±2a, etc.). 

The PS selection technique ensures the very important possibility to discriminate "illegal" 
events: double pulses, a events, etc., and thus to suppress background. An example of a double 
pulse is shown in fig. |]a. Value of the shape indicator for the full signal is SI = -47; for the 
first pulse SIi = 18.4 (hence it corresponds to 7 or (3 particle), for the second pulse SI2 = 37.4 
(a particle). The energy release is 1.97 MeV, and without PS analysis it would be a candidate 
event for 2v2(3 decay of 116 Cd. 

Since the shape indicator characterizes the full signal, it is also useful to examine the pulse 
front edge. For example, it was found that at least 99% of "pure" 7 events (measured with 
calibration 232 Th source) satisfy the following restriction on pulse rise time : At(fis) < 1.24 - 



with conversion electrons of 207 Bi by using the signal of the thin plastic scintillator (placed between 
source and detector) as signature of the electron hitting the 116 CdW04 crystal. 

3 For the 7-s (300-3200 keV) SI 7 = 18.09 - (4.5x 10~ 5 £ 7 ), <j 7 = 2.61 - (4.7x 10~ 4 £ 7 ) + 707/ E 7 , 
while for the a particles (4000-6000 keV) SI a = 29.0; a a = 5.11 - (5.52 xl0~ 4 £ a ) + 5520/£ Q . Here 
all variables are dimensionless (EL and E a are expressed in keV). 
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0.5i? 7 + 0.078-E 2 , where E y is dimensionless variable expressed in MeV. Hence, this filter was 
applied to the background data, and all events, which do not pass the test, were excluded from 
the residual (3/^ spectrum. 

The results of PS analysis of the data are presented in fig. |5|. The initial (without PS 
selection) spectrum of the 116 CdW04 scintillators in the energy region 1.2-4 MeV - collected 
during 4629 h in anticoincidence with active shield - is depicted in fig. |5|a, while the spectrum 
after PS selection of the /3/j events, whose SI lie in the interval SI y — 3.0<r 7 < SI < S7 7 + 2.4<7 7 
and At(/is) < 1.24 - 0.5£ 7 + 0.078^ (98% of (3/^ events), is shown in fig. ||b. From these figures 
the background reduction due to pulse-shape analysis is evident. Further, fig. |c represents the 
difference between spectra in fig. |]a and ||b. These events, at least for the energy above 2 
MeV, can be produced by 228 Th activity from the intrinsic contamination of the 116 CdW04 
crystals (measured by the time-amplitude analysis as described above). Indeed, two decays in 
the fast chain 212 Bi (Q p = 2.25 MeV) -> 212 Po (Q a = 8.95 MeV, T 1/2 = 0.3 fis) -> 208 Pb can 



not be time resolved in the CdWC>4 scintillator (with an exponential decay time ~15 /is [22 
and will result in one event. The example of such an event - recorded by the PS acquisition 
system - is depicted in fig. [|b. To determine the residual activity of 228 Th in the crystals, the 
response function of 116 CdW04 detectors for the 212 Bi — > 212 Po — > 208 Pb chain was simulated 
with the help of GEANT3.21 code and event generator DECAY4. The simulated function 
is shown in fig. ||c, from which one can see that the high energy part of the experimental 
spectrum is well reproduced (x 2 = 1-3) by the expected response for 212 Bi — > 212 Po ^ 208 Pb 
decays^. Corresponding activity of 228 Th inside the 116 CdW04 crystals, deduced from the fit 
in the 1.9-3.7 MeV energy region, is 37(4) /iBq/kg, that is in a good agreement with the value 
determined by the time-amplitude analysis of the chain 220 Rn — > 216 Po — > 212 Pb. Besides, the 
front edge analysis of 80 events with the energy 2.0-4.2 MeV (SI > SI-, + 2.54a 7 ; At > 0.2 /is) 
was fulfilled and the half-life derived from the average time delay between the first and second 



4 The rest of spectrum below 1.9 MeV (fig. ||c) can be explained as high energy tail of the PS selected 
a particles (see fig. ||). 
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part of the signal (see fig. |b) is Ti/ 2 = 0.31(6) /is, in agreement with the 212 Po table value Xi/2 
= 0.299(2) /is ||. 



Fig. |6| represents the spectrum after PS selection of the background events, whose SI lie in 
the interval S7 7 + 2.4<r 7 < SI < SI a + 2.4<j a (~90% of a events). The obtained distribution 
with maximum at 0.95 MeV is well reproduced by the model, which includes all a particles from 
chains in 232 Th and 238 U families. The total a activity of the 116 CdW04 crystals deduced from 
fig. H is 1.4(3) mBq/kg. This value can be adjusted with the activities determined by the time- 
amplitude analysis under usual (for crystals) assumption that secular radioactive equilibriums 
in some chains of 232 Th and 238 U families (like, f. e. 230 Th — > 226 Ra chain) are broken. 

III. RESULTS AND DISCUSSION 

A. Two-neutrino double beta decay of 116 Cd 

To determine the half-life of two-neutrino 2(3 decay of 116 Cd, the background in the energy 
interval 900-2900 keV was simulated with the help of GEANT3.21 package and event generator 
DECAY4. In addition to 116 Cd two neutrino 2(3 decay distribution, only three components 
shown in fig. || were used to build up the background model: 40 K contamination of the enriched 
and natural CdWC>4 scintillators, whose activity limits of less than 4 mBq/kg were established 
earlier [19], and external 7 background caused by 232 Th and 238 U contamination of the PMTs 
(one PMT for 116 CdW04 crystals; one for CdW04; two for plastic active shield)^. This simple 
background model describes experimental data in the chosen energy interval 900-2900 keV 
reasonably well (x 2 =1.3) and gives the following results: the activities of 40 K inside the enriched 
and natural CdWC>4 crystals are equal 0.8(2) and 2.1(3) mBq/kg, respectively; the half-life of 
two-neutrino 2(3 decay of 116 Cd is Ti/ 2 {2i>) = 2.6(1) xlO 19 yr (only statistical uncertainties are 



5 The radioactive impurities of all PMTs used in the installation were previously measured by R&D 
low background set-up as (0.4-2.2) Bq/PMT and (0.1- 0.2) Bq/PMT for 226 Ra and 228 Th activity, 
respectively [plf . 

10 



given, while systematical errors are pointed below). 

Taking advantage of the high statistics in our experiment (approximately 3600 events 
of 116 Cd two neutrino 2(3 decay are contained within the interval 900-2900 keV), we 
can prove our model with the help of experimental 2z/2/3 decay Kurie plot: K{e) = 
[S(e)/ {(e 4 + We 3 + AOe 2 + 60e + SO)^}] 1 / 5 , where S is the number of events with the energy e 
(in electron mass units) in the experimental spectrum after background subtraction. For the 
real 2v2(3 decay events such a Kurie plot should be the straight line K(e) ~ (Q2/3 ~ s), where 
Q 2 p is the 2(5 energy release. From fig. |7|, where our experimental Kurie plot is depicted, one can 
see that in the region 1.1-2.4 MeV it is well fitted by the straight line with Q 2 @ = 2790(87) keV 
(table value is Q 2 p = 2805(4) keV). To take into account the energy resolution of the detector 
the experimental spectrum was also fitted by the convolution of the theoretical 2u2/3 distribu- 
tion p(s) = A ■ e(e i + 10e 3 + 40e 2 + 60e + 30) • (Q 2 p — ^) 5 with the detector resolution function 
(Gaussian with FWHM determined as described above) with A and Q2/3 values as free pa- 
rameters. This fit in the energy region 1.2-2.8 MeV yields a very similar value Q 2 p = 2779(52) 
keV and an A corresponding to Ti/ 2 {2v) = 2.5(3) xlO 19 yr, thus justifying our assumption that 
experimental data in the region above 1.2 MeV are related mainly with 116 Cd two neutrino 2(3 
decay. In fact, a signal to background ratio deduced from our data is 4:1 for the energy interval 
1.2-2.9 MeV, and 15:1 for the energy range 1.6-2.9 MeV, which are higher than those reached 
up-to-date in other 2(5 decay experiments [110. 

To estimate systematical uncertainties of the measured half-life, different origins of errors 
were taken into account, whose contributions are listed in Table |. The final value is equal to: 

T 1/2 (2z/) = 2.6 ± 0.1(stat)±[j;I(syst)xl0 19 yr. 

Our result is in agreement with those measured earlier (2^/2(2//) = 2.6 jlo;§ X 10 19 yr |26 [ 
and Ti/ 2 (2u) = 2.7to;|(stat)to; 9 (syst)xl0 19 yr [14]]) and disagrees to some extent with the value 
T l/2 (2v) = 3.75 ± 0.35(stat)±0.21(syst)xl0 19 yr from ref. f2^Q. 



'Note, that in |27|] the quite small detection efficiency (1.73%) was calculated by the Monte Carlo 
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B. New limits for 0u2/3 decay of b Cd to ground state of b Sn 

To estimate the half-life limit for the neutrinoless decay mode, a simple background model 
was used. In fact, in the 0^2/3 decay energy region only three background contributions are 
important: (i) external 7 background from U/Th contamination of the PMTs; (ii) tail of the 
2v2(3 decay spectrum; and (iii) internal background distribution expected from the 212 Bi — > 
212 Po — > 208 Pb decay ( 228 Th chain). As it was shown above, two decays in the fast chain 212 Bi 
— ► 212 Po — > 208 Pb really create the background in the region of 0z/2/5 decay (see fig. |^c). For the 
activity of 228 Th inside the 116 CdW04 crystals two values were obtained: 38(3) /iBq/kg (time- 
amplitude method) and 37(4) /iBq/kg (pulse-shape analysis). Hence, in the limit of statistical 
errors, we do not find an indication of a failure for the rejection of a pulses by our PS analysis. 

The high energy part of the experimental spectrum of the 116 CdW04 crystals measured 
in anticoincidence with the shielding detectors and after the time-amplitude and pulse-shape 
selection is shown in fig. ||. The peak of 0u2f3 decay is absent, thus from the data we obtain a 
lower limit of the half-life: lim T1/2 = ln2- N-rj-t/ MmS, where N = 4.66 xlO 23 is number of 116 Cd 
nuclei, t is the measuring time (t = 4629 h), r] is the total detection efficiency for Oz/2/3 decay, 
and lim S is the number of events in the peak which can be excluded with a given confidence 
level. The value of the detection efficiency t]mc — 0.83 was calculated by the DECAY4 and 
GEANT3.21 codes, while the efficiency of the PS analysis rjps = 0.98 was determined as 
described above, thus the total efficiency 77 = tjmc • Vps = 0.81. To obtain the value of lim 5, 
the part of the spectrum in the 1.9-3.8 MeV region was fitted by the sum of the simulated 0^2/5 
peak and three background functions (external 7 rays from PMT-s contamination; contribution 
from 212 Bi — > 212 Po — > 208 Pb intrinsic chain; and 2z/2/3 decay tail). This fit gives the value of 
S = — 1.1 ± 1.2 counts, which corresponds - in accordance with Feldman- Cousins procedure 
for results close to the edge of physically accepted area recommended by the Particle Data 



Group (PDG) H] - to a lim S = 0.9(0.3) counts with 90%(68%) C.L., and subsequently to 



method without experimental test, thus perhaps systematical error could be higher than quoted value. 



12 



T 1/2 (0u2(3) > 1.5(4.6) xlO 23 yr at 90%(68%) C.L. However, because of the low statistics in 
the energy range where the effect is expected, the obtained values can be cross-checked in a 
more simple and explicit way. Indeed, in the energy interval 2.6-3.1 MeV (containing 91% 
of 0^2/3 peak) there is only one measured event, while the background expected on the basis 
of the GEANT simulation, in the same energy region is 3.2l 2 ;;| counts (1.9±0.7 events from 
PMT contamination; 0.4±0.1 events from 2v2fi distribution; 0.9lo.g counts from mentioned 
212 Bi — > 212 Po — > 208 Pb chain). Following the PDG recommendation p9|j28[| we can derive 



from these numbers the excluded limit as limS 1 = 1.8(0.5) with 90%(68%) C.L., which leads 
to T 1/2 {0u2(3) > 0.7(2.5) xlO 23 yr at 90%(68%) C.L. confirming the previous estimate. Finally, 
the following values were accepted as conservative half-life limits for neutrinoless 2/3 decay of 
116 Cd: 

T 1/2 (0z/2/3) > 0.7(2.5) xlO 23 yr, 90%(68%) C.L. 



Using calculations |L7|], one can obtain restrictions on the neutrino mass and right-handed 
admixtures in the weak interaction: m v < 3.0 eV, rj < 3.9xl0~ 8 , A < 3.4xl0~ 6 at 90% C.L., 
and neglecting right-handed contribution m u < 2.6(1.4) eV at 90% (68%) C.L. On the basis 



of calculations [^7] we get a similar result: m v < 2.4(1.3) eV at 90%(68%) C.L. In accordance 
with ref. |30| the value of the R-parity violating parameter of minimal SUSY standard model 
is restricted by our T 1/2 limit to e < 8.8(6.4) xl0~ 4 at 90%(68%) C.L. (calculations give 
more stringent restrictions: e < 3.4(2.4)x 10 -4 ). 



C. The bounds on 0v2(3 decay of 6 Cd to excited levels of 6 Sn 

Not only ground state (g.s.) but also excited levels of 116 Sn with Ei ev < Q 2 p can be 
populated in Oz/2/3 decay of 116 Cd. In this case one or several 7 quanta, conversion electrons 
and/or e + e~ pairs will be emitted in a deexcitation process, in addition to two electrons emitted 
in 2(3 decay. The response functions of 116 CdW04 detectors for 0^2/3 decay to the first and 
second excited levels of 116 Sn (2^ with E iev = 1294 keV and 07/ with E [ev = 1757 keV) were 
simulated with the help of GEANT3.21 and DECAY4 codes. The full absorption of all emitted 
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particles should result in the peak with E = Q 2 p (practically the same peak as it is expected for 
the Oz/2/3 decay of 116 Cd to the g.s. of 116 Sn). Calculated full peak efficiencies are: rj(2f) = 0.14 
and rj(0f) = 0.07. These numbers and the value of limS = 1.8(0.5) with 90%(68%) C.L. 
(determined for the g.s.— >g.s. transition) give the following restrictions on half-lives of 116 Cd 
neutrinoless 2(3 decay to excited levels of 116 Sn: 

Ti /2 (g.s.-> 2f) > 1.3(4.8) xlO 22 yr, 90%(68%) C.L, 
r 1/2 (g.s.-> 0t) > 0.7(2.4) xlO 22 yr, 90%(68%) C.L. 

D. Neutrinoless 2/3 decay with Majoron(s) emission 

The procedure to obtain half-life limits for 0v2/3 decay with one (two) Majoron(s) emission 
was carried out in two steps as follows. First, because in the measured spectrum contributions 
of 40 K are negligible above the energy 1.6 MeV, the data were fitted in the energy region 1.6-2.8 
MeV for Oz/Ml mode (1.6-2.6 MeV for 0z/M2) by using only three theoretical distributions: 7 
background from measured PMT-s contamination ( 226 Ra and 232 Th chains) and two neutrino 
2(3 decay of 116 Cd, as background, and Oz/2/3 decay with one (two) Majoron(s) emission, as 
effect. With this simple model the x 2 value was equal to 1.1 both for Oz/Ml and 0z/M2 fits. 
As a result, the number of events under a theoretical Oz/Ml curve was determined as 9 ± 21, 
giving no statistical evidence for the effect. It leads to an upper limit of 41(26) events at 
90%(68%) C.L, that together with an efficiency value rj = 0.905 corresponds to the half-life 
limit 7i/ 2 (0z/Ml) > 3.7(5.9) x 10 21 yr. A similar procedure for 0u2(3 decay with two Majorons 
emission gives T 1/2 (0z/M2) > 5.9(9.4) x 10 20 yr at 90%(68%) C.L. The part of the experimental 
spectrum and theoretical Oz/Ml and 0z/M2 distributions with half-lives equal to these limits 
are shown in fig. [|. On one hand, the obtained results can be treated as conservative because 
the accepted background model consists of only two origins, the external background from 
U/Th contamination of the PMT and 2v2(3 decay distribution, while in the energy region of 
interest some other sources of background, such as the mentioned 212 Bi — > 212 Po — > 208 Pb chain 
contribution from the intrinsic impurities of the crystals, could enlarge the Oz/Ml and 0z/M2 
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limits. On the other hand, the uncertainty of the 2h>2j3 decay half-life of 116 Cd could lead to 
the overestimated value of the OuMl bound. To avoid the last possibility we have estimated 
the Ti/2(0z/Ml) limit in a more explicit way. With this aim the energy interval 2.5-2.8 MeV 
- where the tail of 2z/2/3 decay distribution is practically zero (see fig. ||) - was used as the 
most sensitive region for the OuMl double f3 decay search. In this energy interval the number 
of measured counts is 3, while the expected contribution (in the same energy range) from the 
PMT contamination is 3.2±1.0 events and from 2z/2/3 decay is 1.5±0.5 events, thus the total 
expected background is 4.7±1.1 counts. Following the PDG recommendation [[29| , p8| 



we can 



derive from these numbers the excluded limit for the effect being sought as 3.0 events with 90% 
C.L. Taking into account that the interval 2.5-2.8 MeV contains 8.9% of full Oz/Ml curve, it 
yields a limit of Ti/2 (Oz/Ml) >4.5xl0 21 yr (90% C.L.), confirming the preceding estimate: 

T 1/2 (0vMl) > 3.7(5.9) xlO 21 yr, 90%(68%) C.L., 
T 1/2 (0z/M2) > 5.9(9.4) xlO 20 yr, 90%(68%) C.L. 

Both present half-life limits are more stringent than those established in our previous mea- 



surement during 19986 h |18| and in the NEMO experiment [27]. 

The probability of neutrinoless 2(3 decay with Majoron emission can be expressed as: 
|Ti/ 2 (0z/M1)| l = < g M > 2 • |NME| 2 -G, where < g M > is the effective Majoron-neutrino 
coupling constant, NME is the nuclear matrix element and G is the kinematical factor. Using 
our result T 1/2 (0z/Ml) > 3.7(5.9) xlO 21 yr and values of G and NME calculated in the QRPA 
model with proton-neutron pairing [^| we obtain gM < 12(9.5) xl0~ 5 (gM < 6.5(5.4) x 10 -5 on 
the basis on calculation |2"7] ) with 90%(68%) C.L., which is one of the best restriction up-to-date 
obtained in the direct 2f3 decay experiments 0. 



IV. CONCLUSION 

The search for 116 Cd double (3 decay with enriched 116 CdW04 scintillators has entered in 
a new phase. The set-up with four 116 CdW04 crystals (339 g) is running since October 1998. 
Improved passive shield, new active shield made of fifteen CdW04 crystals (total mass 20.6 
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kg), as well as time-amplitude and pulse-shape analysis of the data result in the reduction 
of the background rate in the 2.5-3.2 MeV region to 0.03 counts/yr-kg-keV. This reduction, 
together with an about threefold increase in the number of 116 Cd nuclei, leads to the substantial 
sensitivity enhancement of the 116 Cd experiment by more than one order of magnitude. Due 
to that the neutrino mass limit of m v < 2.6(1.4) eV at 90%(68%) CL. was set after the first 
4629 h run. 

In August 1999 one of our 116 CdW04 crystals was annealed at high temperature, and its 
light output was increased by ~13%. The PMT of the main 116 CdW04 detectors was changed 
by the special low background EMI tube (5 inches in diameter) with the RbCs photocathode, 
whose spectral response better fits the CdW04 scintillation light. As a result, the spectro- 
metric parameters of four crystals taken as a whole were improved. In particular, the energy 
resolution of the main detector is now 11.4% at 1064 keV and 8.6% at 2615 keV (comparing 
with those before this upgrading: 14.5% and 11%). Besides, the PS discrimination ability of 
the detector was improved too, as it is visible from fig. ^ where the SI distribution of the 
measured background events - before and after the last upgrading - is depicted. It is expected 
that after approximately 5 years of measurements the half-life limit Ti/2(0i/2/3) > 4xl0 23 yr 
will be reached which corresponds to m u < 1.2 eV. The bounds on neutrinoless 2f3 decay with 
Majorons emission and 2(3 transitions to the excited levels of 116 Sn would be improved too. 
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FIG. 1. The scheme of the new set up with four CdW04 detectors. 
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FIG. 2. Background spectrum of 116 CdW04 detectors (339 g) measured in the set up with four 
enriched crystals during 4629 h (solid histogram). The old data obtained with one 116 CdWO"4 crystal 
(121 g; 19986 h) is shown for comparison (thin histogram; the data are normalized to 4629 h and 
mass of the new detector). The background components used for fit in the energy region 900 — 2900 
keV: (a) 2i/2/3 decay of 116 Cd (fit value is T 1/2 (2u) = 2.6(1) xlO 19 yr); (b) 40 K inside the 116 CdW0 4 
detector (activity value from the fit is 0.8(2) mBq/kg); (c) 40 K in the shielding CdWO"4 crystals (fit 
value is 2.1(3) mBq/kg); (d) 226 Ra and 232 Th contamination of PMTs. 
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FIG. 3. The energy spectra of the first (a) and second (b) a particles from the 220 Rn — ► 216 Po 
— ► 212 Pb chain selected by time-amplitude analysis from 116 CdW04 data. Their equivalent energies 
in the /3/'y energy scale are near 5 times smaller because the relative light yield for a particles as 
compared with that for electrons (a//3 ratio) is ~ 0.2. (c) Time distribution between the first and 
second events together with exponential fit (Ti/ 2 = 0.15(1) s, while the table value is T^i-i = 0.145(2) 
s [24]). 
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FIG. 4. (a) Example of a double pulse with total energy 1.97 MeV. The shape indicators for the 
full signal and separately for its first and second parts are Sis = —47; SI\ = 18.4 (close to S7 7 ) and 
SI2 = 37 .4 (close to SI a ). Most probably, this is the couple of successive decays 214 Bi (/?) — > 214 Po 
(a; T 1/2 = 164.3 /js) -> 210 Pb. (b) Probable event of the chain 212 Bi {(3) -> 212 Po (a; T 1/2 = 0.3 ^s) 

_^ 208p] ;) 
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FIG. 5. (a) Initial spectrum of 116 CdWO"4 crystals (339 g, 4629 h) in anticoincidence with shielding 
detectors without pulse-shape discrimination; (b) PS selected P/'y events (see text); (c) the difference 
between spectra in fig. 5a and 5b together with the fit by the response function for 212 Bi — ► 212 Po — ► 
208 Pb decay chain. The fit value is 37(4) /iBq/kg for 228 Th activity inside 116 CdWO"4 crystals. 
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FIG. 6. Spectrum after PS selection of the background events, whose SI lies in the interval 
SIj + 2A<jy < SI < SI a + 2Aa a (it contains ss90% of all a events). The model distribution (smooth 
line) includes all a-particles from chains in 232 Th and 238 U families. The total a activity of the 
116 CdW04 crystals is derived as 1.4(3) mBq/kg. 
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FIG. 7. The 2^2/3 decay Kurie plot and its fit by the straight line in 1100-2400 keV region. 
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FIG. 8. Part of experimental spectrum of the CdW0"4 detectors measured during 4629 h 
(histogram) together with the fit from 2z/2/3 contribution (Tu2 = 2.6xl0 19 y). The smooth curves 
0i/2/3Ml and 0z^2/?M2 are excluded with 90% C.L. distributions of Oz^Ml and 0^M2 decay of 116 Cd 
with Tj/2 = 3.7xl0 21 y and T1/2 = 5.9xl0 20 y, respectively. In the insert the expected peak from 
0u2P decay with r 1 / 2 (0i') = l.OxlO 22 y is shown together with the excluded (90% C.L.) distribution 
(solid histogram) with T 1/2 (0i/) = 7.0 xlO 22 y. 
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FIG. 9. Shape indicator (SI) distributions, which represent the PS discrimination ability of the 
116 CdWO"4 detectors: (a) for background events (0.8-1.2 MeV) collected during 4629 h before the last 
upgrading; (b) for background events (0.8-1.2 MeV) measured during 2734 h after the last upgrading. 
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TABLE I. Different origins of the systematical uncertainties and their contributions to the half-life 
value of 116 Cd two neutrino 2(3 decay 



Origin of the systematical error 




Value range 


Contribution to 
T x / 2 (2v) value, 10 19 yr 


Live measuring time 




96±^ % 


+0.05,-0.2 


Efficiency of PS analysis 




98^ % 


+0.05,-0.3 


Detection efficiency of 2v2f3 decay 
(GEANT model uncertainty) 




96±4% 


±0.1 


90 Sr _90 Y i m p Ur ity in 116 CdW0 4 




<0.17 mBq/kg 


+0.5 


234m Pa impur i ty in H6cdW0 4 




<0.19 mBq/kg 


+0.3 
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